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Abstract: Chiral separations of five [3-adrenergic antagonists (propranolol, esmolol, atenolol, 
metoprolol, and bisoprolol) were studied by capillary electrophoresis using six cyclodextrins 
(CDs) as the chiral selectors. Carboxymethylated-P-cyclodextrin (CM-(3-CD) exhibited a 
higher enantioselectivity power compared to the other tested CDs. The influences of the 
concentration of CM-[3-CD, buffer pH, buffer concentration, temperature, and applied 
voltage were investigated. The good chiral separation of five |3-adrenergic antagonists was 
achieved using 50 mM Tris buffer at pH 4.0 containing 8 mM CM-[3-CD with an applied 
voltage of 24 kV at 20 °C. In order to understand possible chiral recognition mechanisms 
of these racemates with CM-[3-CD, host-guest binding procedures of CM-(3-CD and these 
racemates were studied using the molecular docking software Autodock. The binding free 
energy was calculated using the Autodock semi-empirical binding free energy function. 
The results showed that the phenyl or naphthyl ring inserted in the hydrophobic cavity 
of CM-P-CD and the side chain was found to point out of the cyclodextrin rim. Hydrogen 
bonding between CM-P-CD and these racemates played an important role in the process of 
enantionseparation and a model of the hydrogen bonding interaction positions was 
constructed. The difference in hydrogen bonding formed with the -OH next to the chiral 
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center of the analytes may help to increase chiral discrimination and gave rise to a bigger 
separation factor. In addition, the longer side chain in the hydrophobic phenyl ring of the 
enantiomer was not beneficial for enantioseparation and the chiral selectivity factor was 
found to correspond to the difference in binding free energy. 

Keywords: molecular docking; cyclodextrin; [3-adrenergic antagonists; capillary 
electrophoresis; chiral recognition mechanism 



1. Introduction 

[3-Adrenergic antagonists are effectively used for the treatment of hypertension, prevention of 
anginal attacks, suppression of cardiac arrhythmia, prevention of myocardial infarction and possibly, 
amelioration of congestive heart failure [1]. Most of these drugs are chiral. It is well known that the 
pharmacodynamic and pharmacokinetic properties of these [3-adrenergic antagonist enantiomers differ 
significantly. For example, [3-blocking activity is predominantly related to S-propranolol, which is 
about 100 times more potent than ^-propranolol [2,3]. Therefore, chiral separations of [3-adrenergic 
antagonists are considered to be an essential issue. 

Analytical methods so far used for chiral separations of these [3-adrenergic antagonists include 
thin layer chromatography [4,5], HPLC [6-8], capillary electrochromatography [9,10], and capillary 
electrophoresis (CE) [11-14]. CE especially has attracted the greater interest for chiral separations. It 
has a number of advantages in chiral separations such as high separation efficiency, low solvent and 
selector consumption as well as the ability to readily change the types and concentration 
of chiral selectors [15]. Numerous chiral selectors (e.g., cyclodextrins (CDs), chiral crown ethers, 
proteins, chiral surfactants, macrocyclic antibiotics, ligand-exchange complexes, and polysaccharides) 
are currently available. Among them, CDs are the most widely used due to their excellent chiral 
recognition abilities. Enantioseparation can be achieved using chiral selectors, which discriminate 
between the enantiomers by an enantioselective complexation between the enantiomers of the analyte 
and the chiral selector, giving rise to differences in the electrophoretic mobility of the enantiomers. 
Chiral recognition mechanism, when CDs are used as chiral selectors, is usually based on inclusion 
complexation where the analyte fits into the CD cavity. Although some researchers [16-18] have 
studied chiral recognition mechanisms of [3-adrenergic antagonists with CD by NMR, detailed 
mechanisms underlying the separation remain unknown. The influence of hydrogen bonding in the 
chiral recognition process was especially difficult to explain due to the limitations of the experimental 
methods. Fortunately, molecular modeling methods have been recently proposed as powerful tools to 
obtain information about the emerging interaction of inclusion complexes between CDs and 
enantiomers and then to elucidate chiral recognition processes [19]. In recent years, various molecular 
modeling studies such as PM6 semi-empirical methods [20,21], molecular docking [22-24] and 
molecular dynamics simulations [25-28] have been performed to investigate CD inclusion complexes 
with the aim of comprehending the mechanism of the complexation and to correlate with the 
experimental results. The current study applied Autodock molecular docking technique to gain insight 
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into the selector-enantiomer interaction energy and provided useful supporting information for 
enantiomeric separation. 

In this study, CM-P-CD was used as the chiral selector (target) and five [3-adrenergic antagonists 
(Figure 1) as the ligands. The influences of chiral selector and its concentration, buffer pH, buffer 
concentration, temperature and applied voltage on chiral separation of these enantiomers were studied. 
The course of host-guest inclusion was determined by means of a molecular docking technique and 
then the binding free energy was calculated using the Autodock semi-empirical binding free energy 
function. Based on the simulation data of molecular docking and the experimental results of chiral 
separation by CE, the chiral recognition mechanism of [3-adrenergic antagonists with CM-P-CD was 
studied and the influence of hydrogen bonding in the chiral recognition process was explained. 



Figure 1. The structures of five P-adrenergic antagonists. 




Bisoprolol 



2. Experimental 

2.1. Apparatus 

'IT) 

All CE experiments were performed on an HP CE system from Agilent Technologies (Palo Alto, 
CA, USA) with an online diode-array detector. Instrument control and data acquisition were performed 

'IT) 

by the HP CE ChemStation software [29]. An uncoated fused-silica capillary with a total length of 
48.5 cm (effective length 40 cm) x 50 um i.d. (Yongnian Optical Fiber Plant, Hebei, China) was used 
as the separation column. Before the first use, the new capillary column was rinsed with methanol for 
15 min, followed by 0.1 M hydrochloric acid solution for 15 min, and then activated by flushing with 
1 M NaOH for 30 min. Every day the capillary was rinsed with 0.1 M NaOH for 15 min, water for 
5 min, and conditioned with background electrolyte (BGE) for 15 min. Between each run, it was 
treated with 0.1 M NaOH for 4 min, water for 2 min, and running BGE for 4 min. The samples were 
injected at a pressure of 50 mbar for 3 sec and separated using a constant voltage of 24 kV. The 
detection wavelength was 214 nm and the column temperature was kept at 20 °C. An ORION Model 
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828 pH meter with a precision of 0.01 pH units (Thermo Electron, Milford, MA, USA) was used for 
pH measurement. 

2.2. Chemicals 

Propranolol, esmolol, atenolol, metoprolol, and bisoprolol were purchased from the Chinese Institute 
of Biological Products Control (Beijing, China). P-CD was obtained from Shanghai San Jie 
Biochemistry Technology Co., Ltd. (Shanghai, China). Hydroxypropyl-P-CD (HP-[3-CD), heptakis-(2,6- 
di-0-methyl)-p-CD (DM-p-CD), heptakis-(2,3,6-tri-0-methyl)-p-CD (TM-p-CD) were purchased from 
Sigma-Aldrich. Sulfated-p-CD (S-p-CD) and carboxymethylated-p-cyclodextrin (CM-p-CD) with 
purity higher than 98.0% were synthesized by Department of Organic Chemistry, School of Pharmacy, 
Second Military Medical University (Shanghai, China) and they were identified by IR, NMR, and 
ESI-MS. Tris was purchased from Shanghai Shisheng Cell Biotechnology Co., Ltd. (Shanghai, China). 
Phosphoric acid (H3PO4) was obtained from Shanghai No. 4 Reagent Factory Kunshan Branch 
(Jiangsu, China) and deionized water was prepared by Milli-Q System (Millipore, Bedford, MA, USA). 
All other reagents were of analytical grade. 

2.3. Preparation ofBGE Solutions and Analytes 

The running BGE consisted of 50 mM Tris and 8 mM CM-P-CD unless otherwise stated. The pH of 
the buffer was adjusted with phosphoric acid (H3PO4). 

Stock solutions of five drugs were dissolved in 50 mM Tris buffer at a concentration of 10 mM. The 
samples were prepared by diluting the stock solutions with Tris buffer to a concentration of 0.2 mM 
before injection. 

All solutions were filtered through a membrane filter (0.22 jam) and then stored at 4 °C until analysis. 

2.4. Molecular Construction and Optimization 

Enantiomer construction in its protonated state was carried out using the Sketch Molecule module 
of SYBYL7.0 software package [30]. The molecular mechanics Powell method [31] was applied for 
structure optimization and charges were considered by using the Gasteiger-Huckel charge calculation. 
MOP AC 6.0 [32] was used to adjust the spatial unreasonable bond distances and bond angles of all 
enantiomers [33,34]. 

The three-dimensional structure of CM-P-CD was constructed taking as reference P-CD coordinates 
extracted from the P-CD_a-hemolysin complex crystal structure (PDB entry code: 3M3R [35]). They 
were then optimized using the molecular mechanics Powell method and charges were considered 
by using the Gasteiger-Huckel charge calculation in SYBYL7.0 software package. The structure of 
CM-P-CD was derived from the three-dimensional structure of P-CD. All molecular mechanics 
calculations and quantum chemistry calculations were carried out on an Origin 300 Server. 
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2.5. Docking and Mechanic Calculations 

Autodock 3.0 [36] software was applied for molecular docking and the software had been validated 
for simulating the interactions of two molecules [37]. The optimized drug enantiomers were inserted 
into the hydrophobic cavity of CM-[3-CD. The Lamarckian genetic algorithm in Autodock 3.0 was 
applied to search the conformational and orientational space and the number of iterations for the local 
search was 300 using the Solis and Wets algorithm. The most stable docking conformations were 
selected on the basis of docking energies and the first pose of the most populated cluster. The binding 
free energy AG was calculated from the most stable docking conformations using Autodock 
semi-empirical binding free energy function [38]: 



12 6 

r.. r. 

V v v J 



+ AG hbond Y J E(t) 



r 12 r 10 , 



+ AG tor N tor +AG so! ^(S l V J +S J V i p /2a2 

'J 

where the first three terms are the typical molecular mechanics terms for dispersion/repulsion, hydrogen 
bonding, and electrostatics, respectively; AG tor term models the restriction of internal rotors and global 
rotation and translation; and AG so i term models desolvation upon binding and the hydrophobic effect 
(solvent entropy changes at solute-solvent interfaces). The magnitude of the binding free energy change 
indicates the tendency towards complexation. The more negative the binding free energy change is, the 
more thermodynamically favorable is the inclusion complex. 

2.6. Chromatographic Thermodynamics Calculations 

The equilibrium constant, K\ can be obtained from the binding free energy using Equation (2). 

AG = -RT\nK i (2) 

where R is the gas constant, and T is the absolute temperature. Additionally, the separation factor, 
which is also named selectivity factor, is approximately calculated by the time of the first and second 
migrating enantiomer determined by CE using the Equation (3). 

a = t R lt s ( 3 ) 

Moreover, separation factor can also be defined as Equation (4) [39,40]: 

a = K R IK s ( 4 ) 

Where Kr is the equilibrium constant of i?-enantiomer and K s is the equilibrium constant of 
S-enantiomer, assumed to be K R > K s . Therefore, the absolute value of the difference of the binding 
free energy can be calculated using Equation (5). 

\AAG\ = RTlna (5) 
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3. Results and Discussion 

3.1. Optimization of Chiral Separations 

3 . 1 . 1 . Effect of the Type of CD 

The type of CD plays an important role in enantiomeric separation by cyclodextrin-capillary zone 
electrophoresis. The discrimination power of these six CDs was investigated by means of a preliminary 
screening test. Analysis tests were performed with a 50 mM Tris (pH 3.0) buffer individually 
containing CM-p-CD 8 mM, S-p-CD 8 mM, HP-p-CD 30 mM, DM-p-CD 30 mM, TM-p-CD 30 mM, 
and P-CD 16 mM. Large differences in chiral resolution were obtained for the analytes using these 
different CDs. Separation factors obtained from the chiral separations of five drug racemates with the 
six CDs are presented in Table 1. Only enantiomeric separations of propranolol were observed on 
using P-CD as the chiral selector. Partial enantiomeric separations of several drug racemes were 
achieved with HP-p-CD, DM-p-CD, TM-p-CD and S-p-CD. However, all the drugs could be separated 
using CM-P-CD as chiral selector. CM-P-CD, compared to the other CDs, had the better chiral 
recognition abilities for these enantiomers perhaps due to the carboxymethylated group. Hence, 
CM-P-CD was chosen for the optimization of the separation. 



Table 1. Influence of the cyclodextrin (CD) type on the enantioseparation of five 
P-adrenergic antagonists. 



Enantiomers 






Selectivity Factors (a) 






p-CD 


HP-p-CD 


DM-p-CD 


TM-p-CD 


S-p-CD 


CM-p-CD 


Propranolol 


1.009 


1.016 


1.010 


1.015 


1.019 


1.053 


Esmolol 


1.000 


1.012 


1.007 


1.010 


1.006 


1.036 


Atenolol 


1.000 


1.006 


1.000 


1.000 


1.007 


1.012 


Metoprolol 


1.000 


1.000 


1.000 


1.000 


1.000 


1.010 


Bisoprolol 


1.000 


1.000 


1.000 


1.000 


1.000 


1.007 



HP: hydroxypropyl; DM: heptakis-(2,6-di-0-methyl); TM: heptakis-(2,3,6-tri-0-methyl); 



S: sulfated; CM: carboxymethylated. 

3.1.2. Effect of Buffer pH 

The pH of the buffer is another important parameter affecting chiral separation. Effective charge 
and mobility of the analyte and CM-P-CD directly depends on the pH. Experiments were performed 
using buffer at different pH values in the range of 3.0-6.0 obtained by addition of increasing H3PO4 
amounts to Tris solution. The influence of the pH of the buffer on the enantioseparation of five drugs 
was shown in Figure 2. The variation of the pH had the same effect on the chiral separation factor for 
the five drugs. The chiral separation factor reached a maximum value at pH 4.0. Therefore, the buffer 
of pH 4.0 was selected for the optimal enantioseparation method. 
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Figure 2. Effect of buffer pH on the selectivity factor of five drug enantiomers. Conditions: 
uncoated fused-silica capillary, 48.5 cm (effective length 40 cm) x 50 jam i.d.; 
BGE, 50 mM Tris buffer containing 8 mM CM-p-CD; UV detection at 214 nm; 50 mbar 
for 3 sec; temperature, 20 °C; applied voltage, 24 kV. 
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3.1.3. Effect of the CD Concentration 



The CD concentration is also an essential parameter for the optimization of chiral separations 
because the concentration of the chiral selector affects directly the affinity of the enantiomers for the 
selector [41]. The influence of the concentration of CM-P-CD in the separation buffer on the chiral 
separation of these drugs was investigated in the range from 4 to 10 mM using a 50 mM Tris buffer 
(pH 4.0). As shown in Figure 3, the variation of the CD concentration had also the same effect on the 
chiral resolution for the five drugs. The chiral separation factor reached a maximum value at a 
concentration of 8 mM. A further increase in concentration could result in a slow decrease in the 
separation factor. This trend further confirmed the theoretical model developed by Wren and Rowe [42] 
concerning the existence of a maximum selectivity at a certain concentration of chiral selector. 
Therefore, the 8 mM concentration of CM-P-CD was chosen for further optimization of the separation. 



Int. J. Mol. Sci. 2012, 13 



111 



Figure 3. Effect of the carboxymethylated-P-cyclodextrin (CM-J3-CD) concentration on the 
selectivity factor of five drug enantiomers. Conditions: uncoated fused-silica capillary, 
48.5 cm (effective length 40 cm) x 50 urn i.d.; BGE, 50 mM Tris buffer (pH = 4.0) 
containing CM-[3-CD; UV detection at 214 nm; 50 mbar for 3 sec; temperature, 20 °C; 
applied voltage, 24 kV. 
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3.1.4. Effect of the Buffer Concentration, Voltage, and Temperature 

The effect of the concentration of the buffer on the separation was investigated using Tris buffer 
(pH 4.0) in the concentration range 30-60 mM for the separation. Increase of the buffer concentration 
lowers the electroosmotic flow (EOF) and increases the current and the temperature. With the increase of 
the buffer concentration, the separation factor of five drugs changed slightly. Hence, 50 mM Tris buffer 
was selected for the optimal enantioseparation method. The effect of applied voltage in the range of 
20-30 kV was studied. Application of higher voltage (26-30 kV) resulted in an obvious decrease in the 
separation factor, due to the increased Joule heating. However, use of lower voltage (20-22 kV) resulted 
in longer migration times. So the optimal voltage for analysis was chosen as 24 kV. The effect of the 
temperature in the range of 15-25 °C on the separation factor and migration time of five drugs was 
investigated. The increase in the temperature leads to a reduction of migration times by decreasing the 
buffer viscosity. However, higher temperature could enhance the band broadening and thus decrease 
the peak efficiency and the resolution. Therefore, a temperature of 20 °C was chosen for the 
final method. 



3.1.5. Optimum Conditions for Enantioseparation 



After optimization of the various factors, the final method was as follows: uncoated fused-silica 
capillary, 48.5 cm (effective length 40 cm) x 50 um i.d.; BGE, 50 mM Tris buffer, pH 4.0, with 8 mM 
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CM-P-CD; UV detection at 214 nm; injection, 50 mbar for 3 sec of 0.2 mM solution; temperature, 
20 °C; applied voltage, 24 kV. The separation factors obtained for five P-adrenergic antagonists in 
these optimal conditions were listed in Table 2 and the electropherograms were presented in Figure 4. 

Table 2. Binding free energy of carboxymethylated-P-cyclodextrin (CM-J3-CD) with five 
drug enantiomers and selectivity factors. 



Enantiomers 



Binding Free 
Energy (AG) 
(kcal/mol) 



Difference of Binding Free 
Energy (|AAG|) (kcal/mol) 



Selectivity 
Factors (a) ** 





fl-AG* 


5-AG* 


Calculated 


Experimental 




Propranolol 


-4.340 


-3.926 


0.414 


0.183 


1.078 


Esmolol 


-5.113 


-5.469 


0.356 


0.131 


1.055 


Atenolol 


-5.984 


-5.679 


0.305 


0.091 


1.038 


Metoprolol 


-3.410 


-3.622 


0.212 


0.057 


1.024 


Bisoprolol 


-4.612 


-4.747 


0.135 


0.022 


1.009 



* Average energy of the best cluster (the lowest docked energy); ** Conditions are the same as 
optimum conditions. 

Figure 4. Electropherograms of five P-adrenergic antagonists under the optimal conditions: 
(a) propranolol; (b) esmolol; (c) atenolol; (d) metoprolol; and (e) bisoprolol. 
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3.2. Molecular Docking and Chiral Recognition Mechanism 

Molecular modeling is frequently used to rationalize experimental findings concerning chiral 
recognition by CDs due to the limitations of the experimental methods. The methods are valuable tools 
for obtaining information on the geometry and the interaction energy of the inclusion compounds. 
To understand possible chiral recognition mechanisms of these racemates with CM-P-CD, host-guest 
binding procedures of CM-P-CD and these racemates were studied using the molecular docking 
software Autodock 3.0. The molecular docking configuration of five drugs and CM-P-CD was shown 
in Figure 5. A hydrophobic naphthyl or phenyl ring of these racemates inserted in the hydrophobic 
cavity of CM-P-CD and the side chain were found to point out of the cyclodextrin rim. The molecular 
docking configuration of metoprolol and CM-P-CD had been confirmed by NMR studies [16]. 
Although the electrostatic interaction between the positively charged nitrogens of the side chain and 
the negative charges of CM-P-CD was probably the main chiral recognition interaction [43], the 
intermolecular hydrogen bonding in the chiral recognition process could not be underestimated 
because it could be formed between hydrogen donors or acceptors in the side chain and on the 
cyclodextrin rim. Hydrogen bonding interactions of these racemates with CM-P-CD were summarized 
from the molecular docking studies in Figure 6. Based on the findings from the summary of hydrogen 
bonding and the result of molecular docking studies, a model was generated (see Figure 6). 

Figure 5. The molecular docking configuration between five P-adrenergic antagonist 
enantiomers and CM-P-CD: (a) propranolol; (b) esmolol; (c) atenolol; (d) metoprolol and 
(e) bisoprolol. The left is i?-enantiomer and the right is S-enantiomer. The hydrogen 
bonding is indicated with the green dashed line. 

R-enantiomer S-enantiomer 
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Figure 5. Cont. 




In this model there were three hydrogen bonding interaction position (-NH, -OH and the other 
groups in the side chain) and a hydrophobic interaction position. The hydrophobic interaction between 
the phenyl or naphthyl group and the cavity of cyclodextrin may help in forming stable host-guest 
complexes and the difference in hydrogen bonding may increase chiral discrimination besides 
electrostatic interaction between the positively charged nitrogens of the side chain and the negative 
charges of CM-[3-CD. Particularly, the difference in hydrogen bonding formed with the -OH next to 
the chiral center may help to increase chiral discrimination and give rise to the higher separation factor 
such as propranolol, esmolol, and atenolol. Propranolol was able to be separated using all the CDs. 
Interestingly, propranolol had the highest separation factor of the five drugs regardless of the type 
of cyclodextrin as listed in Table 1 . It might be attributable to the strongly hydrophobic naphthyl group 
of propranolol. The strong plane of naphthalene helped in easily forming a stable inclusion complex. 
Also, naphthalene had a larger structure than the phenyl ring, which gave rise to lower spatial degrees 
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of freedom in the hydrophobic cavity of CM-[3-CD as shown in Figure 5. In addition, the separation 
factors of esmolol, atenolol, metoprolol, and bisoprolol, with a hydrophobic phenyl ring in their 
structures, gradually decreased. Considering the structure of these four drugs, the substitutions on the 
phenyl ring had longer side chains, from esmolol to bisoprolol, which may be unfavorable to host- 
guest binding procedures due to less flexibility and steric hindrance. It is suggested that the longer side 
chain in the hydrophobic phenyl ring of the enantiomer is not good for enantioseparation. All these 
differences in the drug enantiomers probably gives rise to the difference in binding energy in the 
enantiorecognition process and then to the chiral discrimination. 

Figure 6. Model and summary table of hydrogen bonding interaction of 
carboxymethylated-P-cyclodextrin(CM-P-CD) with [3-adrenergic antagonist enantiomers. 

, /CH 3 
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The difference in energies of the inclusion complexes between the enantiomers and CD is probably 
a measure of chiral discrimination, which results in the separation of the enantiomers and the different 
separation factors as observed in the experimental studies. Therefore, the binding free energy AG in 
the course of inclusion between each enantiomer and CM-[3-CD was calculated using Equation (1) and 
the results were listed in Table 2. All CM-P-CD-analyte inclusion complexes had binding free energies 
in the range of -3.410 to -5.984 kcal/mol. The negative values for AG for all the complexes indicated 
the spontaneity of the binding of the guest molecule to the host. The calculated and experimental 
absolute values of the difference of the binding free energy |AAG| were also calculated using the stated 
equation and the results were also listed in Table 2. It showed that the calculated |AAG| magnitudes were 
in the order propranolol > esmolol > atenolol > metoprolol > bisoprolol, for the different structures of 
the side chains. The experimental |AAG| and chiral selectivity factor magnitudes had also the same 
order, which indicated that the chiral selectivity factor corresponded with the difference in binding free 



: R 




HB Position 3 or >= 
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energy. However, the values between the calculated and experimental |AAG| had a significant 
difference in the range of 0.135 to 0.414 and 0.022 to 0.183 kcal/mol, respectively. A number of 
assumptions are usually encountered in modeling chiral separations of chromatographic techniques. 
Such factors as buffer effect, solvation effect as well as entropy difference were usually not considered. 
This could lead to the significant difference between the calculated and experimental |AAG|. In 
addition, the values of experimental |AAG| were also approximately calculated using Equation (5). 
However, the combination of molecular docking and CE experiment results could help us predict the 
enantioseparation of the other structurally related analytes under given experimental conditions. 

4. Conclusions 

Chiral separations of five ^-adrenergic antagonists were achieved by capillary electrophoresis using 
CM-P-CD as the chiral selectors. The good chiral separations of five ^-adrenergic antagonists were 
achieved using 50 mM Tris buffer at pH 4.0 containing 8 mM CM-[3-CD with an applied voltage of 
24 kV at 20 °C on investigating the influences of several important analytical conditions. The 
host-guest binding procedures of CM-[3-CD and these drug racemates were studied and their binding 
free energies were calculated using the molecular docking software Autodock 3.0. According to the 
simulation data of molecular docking combined with the experimental results of chiral separation by 
CE, the results achieved were explained as follows: (i) Hydrogen bondings between CM-P-CD and 
these racemates play an important role in the process of enantionseparation; (ii) A longer side chain in 
the hydrophobic phenyl ring of enantiomer was not good for enantioseparation; (iii) The chiral 
selectivity factor corresponded to the difference in binding free energy. It should be noted that a 
different enantio-recognition mechanism may be achieved using the other CDs as the chiral selector. 
These findings are only related to CM-P-CD as the chiral selector. Although it cannot completely 
elucidate chiral recognition mechanisms based alone on CE separation data and molecular modeling 
results, the molecular modeling technique, however, provides us with a good perspective of 
enantio-separation and serves as a useful method for studying chiral recognition mechanisms. 
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